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ON 




(57) Abstract: The invention concerns the 
modelling of interactions between an incident 
wave and an obstacle, in particular in the field 
of non-destructive control. The invention is 
characterized in that it consists in: meshing 
the obstacle surface and assigning at least 
one source (S-,) to each surface element 
(dSi); then calculating conditions at the 
boundaries in each mesh of the obstacle and 
deducing therefrom source values; based on 
an interaction matrix and said source values, 
estimating a physical quantity representing the 
interaction between the wave and the obstacle, 
in any point of space. 

(57) Abrege : U invention conceme la mod61i- 
sation des interactions entre une onde incidente 
et un obstacle, notamment dans le domaine du 
contrdle non destructif. Selon Tinvention, on 
maille la surface de Tobstacle et Ton attribue 



au moins une source (Si) ^ chaque 616ment de surface (dSi). On calcule ensuite des conditions aux limites en chaque maille de 
Tobstacle et Ton en d^duit des valeurs de source. A partir d'une matrice d'interaction et de ces valeurs de sources, on estime une 
grandeur physique representative de Tinteraction entre Tonde et Tobstacle, en tout point de Tespace. 
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Procgde pour gvaluer une grandeur physique reprgeentative 
d'une interaction entre une onde et un obstacle 

L' invention concerne la modilisation des interactions 
5 entre une onde incidente et un obstacle de cette onde, 
notatnment dans le domaine du controle non destructif- 

La mod^lisation des interactions entre une onde et un 
obstacle recevant cette onde, tel qu'une cible plac^e dans 
10 la zone sensible d'un capteur, trouve une application 
avantageuse en contr61e non destmctif . 

On connalt un proc6d6 de modSlisation dit "par ^l^ments 
finis" consistant a appliquer un pavage de I'espace 
15 tridimensionnel entourant 1' obstacle et k ^valuer les 
interactions pr^citSes pour tous les pavSs de I'espace. 

Les proc§d6s de calcul par "€l§ments finis" apportent une 
resolution > un probletne pos6 sous forme d' 6c[uations 
diff^rentielles partielles. lis se fondent sur une 
representation de I'espace d' etude par un assemblage 
d' elements finis, t I'interieur desquels des fonctions 
d' approximations sont determinees en termes de valeurs 
nodales de la grandeur physique recherchee. Le probldme 
physique continu devient done un problSme discret aux 
elements finis oiX les grandeurs nodales sont les nouvelles 
inconnues. De tels precedes cherchent done ^ approcher la 
solution globale, plut6t que les equations axix derivees 
spatiales partielles de depart. 



20 
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30 
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La discretisation de I'espace pris en compte assure que ce 
dernier soit entidrement recouvert par des elements finis 
(lignes, surfaces nn volumes), cette operation porte le 
nom de "maillage" dans I'espace bidimensionnel (2D) ou de 
5 "pavage" dans I'espace tridimensionnel (3D). Les elements 
mis en jeu sont soit rectangulaires ou triangulaires en 
2D, soit paraliei^pip^diques ou t^traedriques en 3D. lis 
peuvent itre de tailles differentes, r^partis \anif ormement 
sur la surface ou non. 

10 

En general, la grandeur physique recherch^e, telle qu'xan 
potentiel electrostatique ou une valeur de pression, est 
connue sur la frontiSre du domaine. Cette frontidre peut 
Stre fictive. On y impose des conditions aux limites. Le 
15 potentiel est done inconnu k I'int^rieur du m§me domaine. 
On dSfinit alors un nceud comme 6tant un sommet d'un 
Pigment. Les inconnues du probldme sont done les valeurs 
du potentiel en chaque noeud de 1' ensemble du domaine. 

20 A titre d' illustration, la figure 6 de I'art ant^rieur 
repr^sente un exemple de surface, constitute par deiix 
mattriaux Ml et M2, de proprittes electromagnttiques 
differentes, et maillee par des Elements triangulaires 
comportant chacun trois noeuds Ai, Bi et Ci. L' ensemble du 

25 domaine est delimit^ par une front iere F. 

Une fois le maillage defini, plusieurs approches existent 
pour transformer la formulation physique du probldme en 
une modeiisation discrete aux Elements finis. Si le 
30 probldme est formulS par des Equations dif f Srentielles et 
consiste k minimiser une f onctionnelle, on applique 
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genferalement une methode variationnelle. Cette 
transformation aboutit k une formulation matricielle dont 
la resolution donne les solutions nodales, les solutions 
aux points non nodaxax 6tant obtenues par interpolation 
5 linSaire. 

Neanmoins, de tels calculs, en trois dimensions, 
requiferent des ressources inf ormatiques importantes et 
ggnerent des temps de calcul tres longs, en depit de 
10 1 * accroissement des performances de logiciels permettant 
d' impl^menter ces calculs. 

Certes, les problSmes 2D, souvent simplifies par des 
conditions de symStrie avantageuses pour ne modSliser 

15 qu'une partie de la gSometrie, se resolvent rapidement. 
Mais il n'en est rien pour les probldmes 3D, qui sont les 
plus frequents. La figure 6 montre combien la finesse du 
maillage, c'est & dire le rapport entre la taille d'un 
element et celle du plus petit detail du domaine, fait 

20 cro£tre le nombre de noeuds. 

En consequence, la quantity d' equations et d'inconnues 
augmente proportionnellement , et, de 1^, le temps de 
calcul n^cessaire k la resolution du probleme. Il est 
25 important de preciser que la generation du maillage, a 
savoir la discretisation de I'espace de travail, et la 
generation de la liste des noeuds consomme un temps de 
calcul superieur k celui necessaire S la resolution du 
probldme . 

30 

La presente invention vient ameiiorer la situation • 
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4 

Elle propose a cet effet tin proc6d6 pour ^valuer vine 
grandeur physique associSe a une interaction entre une 
onde et uii obstacle, dans une region de I'espace 
5 tridimenslonnel , dans lequel : 

a) on determine, par maillage, une plurality 
d'#chantillons de surface dont une partie au moins 
repr€sente la surface d'un obstacle recevant une onde 
principale et 6mettant, en rSponse, une onde secondaire, 
et I'on attribue k chaque ichantillon de surface au moins 
une source imettant \ane onde §16mentaire repr^sentant line 
contribution a ladite onde secondaire, 

b) on forme un systdme matriciel comport ant : 

- une matrice d' interaction, inversible, appliquSe a 
ime region donnSe de I'espace et comportant \in nombre 
de colonnes correspondant a un nombre total de 
sources , 

- une premiere matrice colonne dont chaque coefficient 
est associg k une source et caract^rise I'onde 
#14mentaire qu'elle €met, 

- et une seconde matrice colonne, obtenue par une 
multiplication de la premiere matrice colonne par la 
matrice d' interaction et dont les coefficients sont 
des valeurs d'lme grandeur physique representative de 
I'onde gmise par 1« ensemble des sources en ladite 
region donn^e, 

c) pour estimer les coefficients de la premiere matrice 
colonne, on affecte des valeurs de grandeur physique 
choisies a des points pr6d6termin6s , associ^s chacun a lui 
€chantillon de surface, lesdites valeurs choisies 6tant 
rang€es dans la seconde matrice colonne, et I'on multiplie 
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cette seconde matrice colonne par 1" inverse de la matrice 
d' interaction appliquSe aixxdits points pr€d6tenninSs, 
d) pour ^valuer la'^ite grandeur physique reprisentant 
I'onde 6mise par 1" ensemble des sources dans une region 
5 donn€e de I'espace tridimensionnel, on applique la matrice 
d' interaction ^ ladite region donnee et on multiplie cette 
matrice d' interaction par la premiere matrice colonne 
comportant les coefficients estimes a l'6tape c) • 

10 Ainsi, selon I'un des avantages que procure la presente 
invention, l»§tape de maillage a) ne conceme qu'une ou 
plusieurs surfaces, tandis que le precede de mod^lisation 
de type "par elements finis" n§cessite un pavage de tout 
I'espace avoisinant 1 'obstacle, ce qui permet de reduire, 

15 dans la mise en oeuvre du proc6d6 selon 1' invention, les 
ressources mgmoires et les temps de calcul nfecessaires. 

Le procSd§ selon 1' invention s' applique aussi bien k une 
onde principale Smise par une source lointaine qu'S une 
20 onde principale Smise en champ proche, 

Avantageusement , pour evaluer une grandeur . physique 
representative d'une interaction entre un element 
rayonnant une onde principale en champ proche et un 

25 obstacle recevant cette onde principale, 

- a I'etape a), on determine en outre, par maillage, une 
plurality d' echanti lions de surface reprSsentant ensemble 
une surface active de l'616ment rayonnant I'onde 
principale, et I'on attribue k chaque #chantillon de la 

30 surface active au moins une source ^mettant une onde 
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Sl^mentaire representant une contribution k ladite onde 
principals, 

- on applique les ^uapes b) , c) et d) avx echantillons de 
la surface active, et 

5 - on lvalue ladite grandeur physique representant 
1 ' interaction entre 1* element rayonnant et 1' obstacle dans 
une region donnee de I'espace tridimensionnel , en prenant 
en compte la contribution, en ladite region donnee, de 
I'onde principals 6mise par 1' ensemble des sources de la 

10 surface active et la contribution de 1 ' onde secondaire 
§mise par 1' ensemble des sources de la surface de 
1 'obstacle. 

On entend par les termes "Sl^ment rayonnant" aussi bien un 
15 §metteur de I'onde principale, tel qu'un gSnerateur 
d'ondes, qu'tin r§cepteur de I'onde principale, tel qu'im 
capteur de cette onde, 

Dans un premier mode de realisation, la grandeur physique 
20 ^ evaluer est une grandeur scalaire et, a 1 » 6tape a), on 
attribue une source -unique a chaque gchantillon de 
surface, 

Dans un second mode de realisation, la grandjeur physique a 
25 €valuer est Une grandeur vectorielle exprimSe par ses 
trois coordonn€es dans I'espace tridimensionnel, et on 
attribue, S. l'6tape a), trois sources a chaque echantillon 
de surface. 

30 Dans une realisation avantageuse, pour estimer, Sl I'Stape 
d) , la contribution de I'onde secondaire dans la region 
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doimSe de I'espace, les valeurs de grandeur physique 
choisies k l"§tape c) sont fohction d'un coefficient 
predetermine de reflexion et/ou de transmission de I'onde 
principale par chaque #chantillon de surface de 
5 1' obstacle. 

Ainsi, on comprendra que I'onde secondaire peut aussi bien 
correspondre a une reflexion de I'onde principale, qu'a 
une transmission de l*onde principale, ou encore S une 
10 diffraction de I'onde principale. Dans cette realisation 
avantageuse, I'^tape c) correspond finalement ^ une 
deteinnination des conditions aux limites k la surface de 
1' obstacle, en tant cju ' interface entre deux milieux 
distincts notamment dans une h6t§rostructure . 

15 

Par ailleurs, pour un contr61e non destructif d'une cible 
formant un obstacle d'une onde principale, on attribue un 
coefficient de reflexion ou de transmission choisi k tous 
les points predetermines de la surface de la cible, et on 

20, compare une simulation obtenue par la mise en oeuvre du 
precede au sens de 1' invention avec une mesure 
experiment ale. Ainsi, les points de la surface de la cible 
qui, dans la mesure experimentale, ne verifient pas la 
simulation correspondent a des inhomogeneites ou a des 

25 impuretes sur la surface de la cible. 

Dans une autre approche, on connalt les proprietes 
globales de 1' obstacle, notamment en transmission et/ou en 
reflexion. Par la mise en ceuvre du precede de 1» invention, 
30 on optimise alors la position dans I'espace d'un capteur 
ou meme la forme de ce capteur, pour une application k un 
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contr61e non destructif, ce capteur Stant destinS a 
analyser xrne cible formant un obstacle de I'onde 
principale . 

A cet effet, dans tone realisation avantageuse, on compare 
une plurality de valeurs de la grandeur physique estim^e a 
I'Stape d) du proc#d6 au sens de 1' invention, obtenues 
pour une plurality de regions de I'espace, pour 
s61ectionner une region candidate pour la disposition d'un 
Pigment rayonnant desting a interagir avec 1' obstacle. 

Cotnme indique ci-avant, on entend par les termes "Pigment 
rayonnant" aussi bien un capteur qu'xm ggnSrateur de 
I'onde. On comprendra ainsi que 1 ' optimisation de la 
position de 1' Pigment rayonnant peut s'appliquer aussi S 
1 'optimisation de la disposition ou de la forme d'un 
g^ngrateur d'ondes. Par exemple, la prgsente invention 
trouve encore une application avantageuse & la disposition 
de haut-parleurs dans im volume fermg, d61imit6 par des 
obstacles, comme par exemple I'habitacle d'un v6hicule 
automobile . 

D'autres caractSristiques et avantages de 1' invention 
apparaitront k I'examen de la description d#taill6e ci- 
aprSs, et des dessins annexes sur lesquels : 
- la figure lA repr^sente sch^matiquement les surfaces 
respectives d'un Pigment rayonnant ER gmettant une onde 
et d'un obstacle OBS recevant cette onde, maillges en 
vue d'ivaluer une grandeur scalaire representative de 
I'onde en un point M de I'espace tridimensionnel ; 
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- la figure IB reprSsente en detail iin #chantillon de 
surface dSi correspondant a une* maille de la figure lA, 
ainsi qu'une source Si associSe a I'echantillon de 
surface dSi ; 

5 - la figure 2A repr€sente schematiquement les surfaces 
respect ives d'lin element rayonnant ER gmettant une onde 
et d'un obstacle OBS recevant cette onde, maillees en 
vue d'^valuer une grandeur vectorielle representative 
de I'onde en un point M de l"espace tridimensionnel ; 

10 - la figure 2B reprgsente en detail un echantillon de 
surface dS± correspondant a une maille de la figure 2A, 
ainsi que trois sources associSes SAi, SBi et SCi ; 

- la figure 2C repr§sente, en vue de face, une surface 
maill^e dont chaque Echantillon de surface comporte 

15 trois sources SAi, SBi et SCi, pour 1 'estimation d'une 

grandeur vectorielle ; 

- la figure 3A represente, a titre illustratif, les 
armatures d'un condensateur , de potentiels electriques 
respectifs VI et V2, pour 1* estimation d»un potentiel 

20 glectrique au point M de I'espace tridimensionnel, a 

chaque echantillon de surface dSi de la figure 3A 6tant 
associee une unique source Si ; 

- la figure 3B reprSsente, a titre illustratif, les 
armatures d'un condensateur, de potentiels Electriques 

25 respectifs VI et V2, pour 1 'estimation d'un champ 

Electrique E(M) , au point M de I'espace 
tridimensionnel, a chaque Echantillon de surface dSi de 
la figure 3B Etant associEe trois sources SAi, SBi et 
SCi ; 
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- la figure 4A repr^sente, comme les figures lA et 2A, 
une interaction entre un Pigment rayonnant ER et un 
obstacle OBS, pour 4.r«iuer une grandeur physique 
(scalaire ou vectorielle) en un point M dans une 
5 portion de I'espace dSlimitSe par la surface de 

I'^l^ment rayonnant et la surface de 1' obstacle, ce 
point de I'espace M recevant a la fois I'onde 6mise par 
1- Pigment rayonnant et I'onde rgflichie par 
1 ' obstacle ; 

10 - la figure 4B, complement aire de la figure 4A, 
reprgsente une transmission par 1- obstacle OBS de 
I'onde gmise par 1- Pigment rayonnant ER, en un point M 
d'un demi-espace delimits par le plan que forme la 
surface de 1' obstacle OBS ; 

- la figure 5A reprgsente sch^matiquement un obstacle 
OBS, de dimensions finies, avec des sources associ6es 
aux Schantillons de surface agencies pour estimer une 
grandeur representative d'une rSflexion de I'onde sur 
1 ' obstacle ; 

- la figure 5B, en complement de la figure 5A, represente 
schgmatiquement un obstacle OBS, de dimensions finies, 
avec les sources associees aux 6chantillons de surface 
disposees pour estimer une grandeur representative de 
la transmission de I'onde par 1 'obstacle ; 

- la figure 5C reprgsente une simulation d'une onde 
ultrasonore gmise par un element rayonnant ER et se 
propageant vers un obstacle OBS ; 

- la figure 6 reprSsente un maillage de milieux 
tridimensionnels, pour 1 'application d'un precede de 
calcul par "SlSments finis", au sens de I'etat de la 
technique ; 



20 



25 
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la figure 7A repr^sente en detail vtn Pigment de surface 
et \in point d' observation M dont les positions 
relative pont repSrSes par un an§le G ; et 
la figure 7B repr€sente schematiquement une surface a 
mailler de forme complexe, en particulier avec un point 
d' observation M se situant dans une zone d' ombre par 
rapport a certaines sources de la surface. 



On se rgfdre tout d'abord & la figure lA, sur laquelle la 
10 surface d'un obstacle OBS, recevant xine onde, est maill^e 
selon une plurality d' gchanti lions de surface dSi k dS*, 
conformSment a l'€tape a) prScitge. 

En se rgfgrant ^ la figure IB, a chaque 6chantillon de 
15 surface dSi est associg un h#misphdre HEMi, tangent t 
1 ' gchantillon de surface dSi en un point de contact Pi. 
PrgfSrentiellement, ce point de contact Pi correspond au 
sommet de I'h^misph^re HEMi. Pour 1' estimation d'une 
grandeur physique scalaire au point M (tels qu'xin 
20 potentiel ilectrostatique, une pression acoustique ou 
autre), une source unique Si est associee I 1 ■ ichantillon 
de surface dSi. Comme on le verra plus loin, dans le cas 
de 1' estimation d'une grandeur vectorielle dans un point 
de I'espace M, on affectera plut6t trois sources k chaque 
25 ^chantillon de surface dSi. 



Prgfgrentiellement, I'hgmisphgre HEMi est construit comme 
d^crit ci-aprds. Pendant I'^tape de maillage a) pr6cit6e, 
on lvalue, d'une part, la surface de 1' obstacle OBS et, 
d' autre part, on choisit un nombre d' gchanti lions de 
surface dSi selon la precision souhaitSe de 1' estimation 
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de la grandeur physique au point M. Ainsi, la surface d'un 
gchantillon dSi est donnSe par So/N, oH So correspond k la 
surface totale de 1' obstacle et N correspond au notnbre 
choisi d"€chantillons de surface dSi, 

5 

L»h6misphdre HEMi est de mSme surface que 1 ' Schantillon 
dSi. Ainsi, le rayon Ri de 1 » hemisphere se deduit de 
.1 ' expression : 

^ N 

10 

Chaque maille que represente un ^chantillon de surface dSi 
prSsente, dans I'exemple d^crit, une forme de 
parall61ograihme, de centre Pi correspondant au point 
d' intersection des diagonales de ce parallSlogramme . 
15 L'hgmisphgre HEMi est tangent k 1 ' echantillon de surface 
dSi en ce point Pi. Bien entendu, les mailles peuvent §tre 
de forme diffSrente, triangulaire ou autre. On indique de 
fa?on g€n§rale que le point Pi correspond au barycentre de 
la maille. 

20 

On dgfinit ainsi la position de la source Si (situ^e au 
centre de I'h^misphdre HEMi). La distance sgparant la 
source Si du point de contact Pi correspond au rayon Ri de 
I'h^misphdre HEMi et la droite qui passe par les points Pi 
25 et Si est orthogonale k la maille dSi. 



Dans I'exemple represente sur la figure lA, on maille en 
outre la surface d'un Element rayonnant ER, correspondant 
par exemple k un gSnSrateur d'ondes. A chaque maille de la 
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surface de 1' element rayonnant est associg \m 6chantillon 
de surface dS'i, comme on le verra plus loin. 

Le syst^me matriciel que I'on met en forme k l'6tape b) 
5 precitSe correspond a : 



^V(Mi) ^ 




hi] 


V(M2) 


= F X 


V2 


.V(Mn); 







[1] 



ou 



10 



15 



20 



25 



les coefficients Vj (avec j=l,2,...,N) de la premidre 
matrice colonne correspondent a des valeurs associSes 
chaciine a iine source, telle qu'une charge glectrigue 
(dans le cas de 1 ' estimation d'un potentiel 
#lectrique) , ou encore k un flux magnStique (dans le 
cas de 1' estimation d'un potentiel magnet ique ) , ou 
encore a une vitesse sonore (dans le cas de 
1' estimation d'une pression acoustique liSe k la 
propagation d'xine onde sonore) ; 

les coefficients V(Mi) (avec i=l,2,,..,N) de la seconde 
matrice colonne correspondent chacun k une valeur de la 
grandeur physique (un potentiel §lectrique ou 
magnetique ou une pression) k estimer en un point Mi de 
1 ' espace ; 

la matrice d' interaction P comporte des coefficients 
Ci,j dont 1 'eajpression gSnSrale est donn€e par : 



(Xj=f{MtSj) 



[2] 
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On comprendra ainsi que les coefficients de la matrice F 
sont des coefficients d ' interaction qui dependent de la 
distance s6parant chaque point de I'espace Mi d»une source 
Sj associSe k une thaille dS j . 

Dans le cas de la propagation d'une onde glectrique, les 
coefficients Ci,j, Vj et V(Mi) , respectivement de la 
matrice d' interaction de la premiere et de la seconde 
matrice colonne, sont donnes par : 

^ - 1 

27^Q MiSj 

Vj =qj [3] 

V(Mi) = Ui 

oil : 

- So est une constante dlSlectrigue, 



MjSj est une distance mesuree en valeur alg^brique, 

15 - <Ij correspond k une charge Slectrique caract^risant lane 
source Sj, et 

- Ui correspond S un potentiel 61ectrique au point Mi. 

Dans le cas de la propagation d'une onde magnStique, 
20 1' expression de ces coefficients est la suivante : 



r - ^ 



2jrMo MiSj 

Vj =<Pi C4] 



V(Mi) = 0i 
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Oil : 

- fiQ correspond S la pexm^abiiice magnStique du milieu oil 
se situe le point Mi, 
5 - (pj correspond au flux magn^tique associ^ a la source 
Sj ; 

9i correspond au potentiel magnSticjue au point Mi. 

Dans le cadre de la propagation d'une onde ultrasonore, 
10 ces coefficients sont donnSs par : 

^ icop exp(ikM^) 



MiSj 



Vj=Vj 



[5] 



V(Mi) = pi 

dans laquelle : 
15 - i2="l. 

CO est la pulsation de I'onde sonore ; 

- p est la density du milieu dans lequel se situe le 
point Mi ; 

le vecteur vj correspond & la vitesse sonore issue de 
20 la source Sj ; 

k correspond au vecteur d'onde de I'onde sonore ; et 

- Pi correspond & la pression acoustique gSn^rSe par la 
propagation de I'onde ultrasonore au point Mi. 
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Dans 1' expression des coefficients Ci.j, le terroe dSj 
correspond 4 la surface de 1 '^chantillon associ6 §l la 
source Sj . De pr^f^rence, le maillage d'une surface au 
sens de I'gtape a) du proc6d4 selon 1' invention, est 
choisi de sorte que chaque maille comporte une m§me 
surface dS=dSi=dS2=...=dSj . 

On remarque en particulier dans 1' expression des 
coefficients ci,j qu'ils dependent du produit scalaire 
entre le vecteur d'onde fet le vecteur Mjs] . Ainsi, pour 
des ondes ultrasonores , on tient compte d'un dgphasage 
entre les chemins qui lient chaque source Sj ^ un point de 
I'espace tridimensionnel M, ce dSphasage gtant dfi ^ vine 
difference de marche entre les rayons partant de chaque 
source et arrivant au point M (comme le montre la figure 
4A) . En particulier, on comprendra que 1' angle d' incidence 
d'un tel rayon est pris en compte dans 1 'expression des 
coefficients de la matrice d' interaction F. 

Bien entendu, dans le cadre de la propagation d'une onde 
eiectromagn^tique de frequence 61ev€e, done de longueur 
d'onde courte, qui diffdre du cadre glectrostatique ou 
magnetos tatique ci-avant, on tient compte du terme de 
propagation exp^£r) dans 1' expression de la matrice 
d' interaction, par rapport a la ggom^trie du probldme a 
rSsoudre, comme dans le cadre de la propagation d'vme onde 
ultrasonore ci-avant (relation [5] ) . 

Ainsi, le syst^me matriciel de 1' Equation [1] permet 
d'estimer, S partir d'une matrice d' interaction F et d'un 



wo 2004/044790 



PCT/FR2003/003323 



17 

vecteur comportant les valeurs Vj associSes aux sources 
Sj, les coefficients d'un vecteur (matrice colonne) 
comportant les valeurs de grandeur physique V(Mi) aux 
points de I'espace Mi. 

Pour determiner les valeurs des sources Vj, on applique le 
systdme matriciel suivant : 







fvi \ 


V(P2) 


= F X 


V2 


.V(Pn), 







[6] 



10 

ou : 

- les coefficients de la matrice d' interaction F 
s'expriment par avec Gj=f{liSj) 

- les indices i et j correspondent respect ivement a la 
15 i^"*® ligne et la j^™® colonne de la matrice d' interaction 

P. Cette matrice d» interaction comporte, pour la 
dStemnaination des valeurs associees aux sources Vj, N 
lignes et N colonnes, en rappelant que N est le nombre 
total de mailles k la surface de 1* obstacle ; 
20 - les points Pi correspondent au sommet des hemispheres 
HEMi de la figure IB. 

La mise en oeuvre de I'Stape c) du proc€d6 au sens de la 
prgsente invention correspond a calculer une condition aux 
25 limites pour les points Pi, de propri€t6s connues, comme 
on le verra plus loin. 
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Le systetne matriciel de 1" Equation [6] devient alors : 



fV(Pi) ^ 
V(P2) 

IV(Pn)J 



[7] 



5 oH : 

- F"^ correspond k 1' inverse de la matrice d' interaction 
F ; et 

les valeurs V(Pi) sont predeterminees , en fonction des 
conditions aux limites pr6cit6es. 

10 

On determine ainsi les valeurs de source Vj . 

A partir de 1 " estimation de ces valeurs de source Vj, on 
peut calculer la grandeur physique scalaire en un point M 
15 quelconque de l*espace tridimensionnel, a partir de la 
relation : 



W==ly(WS>>y [8] 

20 Pour obtenir cette expression de la grandeur scalaire 
V{M), la matrice d* interaction F peut ne comporter qu'une 
ligne de coefficients Cj, avec : 

Cj = f (MSj) , 
mais comporte tou jours N colonnes. 

25 

En se r^fSrant Sl nouveau §l la figure lA, on comprendra que 
la surface de 1" obstacle OBS recevant I'onde qu»§met 
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1' Pigment rayoimant ER agit, elle-m§me, comme une surface 
active rg-gmettant une onde secondaire (par example par 
reflexion) . Chaque source Si repr^sente une contribution a 
1' Emission de cette onde secondaire. 

5 

En outre, pour tenir compte k la fois de la presence de 
I'onde principale . et de la presence de 1-onde secondaire 
au point M, on estime la contribution de I'onde principale 
et la contribution de I'onde secondaire au point M par le 
10 systSme matriciel : 





fvi ^ 






V(M) = Fx 


V2 


+ Fx 


V2 

• 











[9] 



oil ; 

- P' est la matrice d' interaction entre la surface de 
I'SlSment rayonnant et le point M ; 

- Vj (avec j=l,2,3,...,N') est la valeur des sources 
attributes k chaque tchantillon de surface dS'j de 
1* fitment rayonnant, N' #tant le nombre total de 
mailles choisi pour la surface active de 1« element 
rayonnant ER. 



Les coefficients de la matrice P' sont encore fonction de 
la distance MS'j, oil S'j sont les sources affecttes k 
chaque tchantillon de surface dS'j de 1' Pigment rayonnant. 

Selon une caractSristique avantageuse, les valeurs des 
sources de 1' obstacle Vj sont d6termin6es en fonction des 
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valeurs des sources de l'gl6ment rayonnant v'j, lesquelles 
sont elles-m§mes calcul^es comme on le verra plus loin en 
reference aux figures 4A, 4B, 5A et 5B. 

5 On se rgfdre maintenant a la figure 2A, dans laquelle on 
affecte trois sources k chaque 6chantillon de surface dSi, 
en vue d'estimer une grandeur physique vectorielle V(M) , 
en un point M de 1 ' espace tridimensionnel . 

10 En effet, on comprendra que pour estimer la grandeur 
vectorielle, par ses trois coordonnSes dans 1' espace x, y 
et z, le nonODre d' Equations a resoudre par rapport au 
systdme matriciel pr4c6dent doit gtre multiplig par trois, 
Ainsi, la matrice F'^ de la relation [7] doit comporter 

15 trois fois plus de lignes que precedemment . La matrice 
d' interaction F doit, elle-meme, comporter trois fois plus 
de colonnes que precedemment et, a cet effet, on pr^voit 
avantageusement trois sources par maille lorsqu'il s'agit 
de determiner les coordonnSes dans 1' espace 

20 tridimensionnel d'un vecteur V(M) . 

En se referant a la figure 2B, les trois sources SAi, SBi, 
SCi, attributes a un tchantillon de surface dSi sont de 
positions respect ives dttermintes comme on I'indique ci- 

25 aprds. Telles que representees sur la figure 2B, les trois 
sources SAi, SBi, SCi sont coplanaires et le plan 
comportant ces trois sources comporte en outre la base de 
1» hemisphere HEMi. L' hemisphere HEMi est construit comme 
indiqu6 ci-avant (de meme surface que la surface de la 

30 maille), avec toutefois le centre de 1« hemisphere qui 
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correspond ici au barycentre des trois sources SAi, SBi et 
SCi. 

On entend par "centre de 1 "hemisphere" le centre du disque 
5 qui constitue la base de 1 •hemisphere . 

Les trois sources qui sont attributes a 1 » tchantillon de 
surface dSi sont disposees aux sommets d'un triangle 
equilateral dont le barycentre Gi correspond au centre de 
1 'hemisphere. De preference, chaque source SAi, SBi et SCi 
est disposee sur le milieu d'un rayon Ri de 1 'hemisphere. 
Ainsi, ies droites qui relient le barycentre Gi a chaque 
source sont ecartees angulairement de 120®. 

En se referant a la figure 2C, 1 » orientation angulaire des 
triangles que ferment les triplets de sources est choisie 
aieatoirement, d'un echantillon de surface t 1» autre. 
Avantageusement , on evite ainsi des artefacts de sur- 
periodicite, dans 1' estimation de la grandeur vectorielle 
au point M, qui pourraient resulter du choix d'une m§me 
orientation angulaire de ces triangles. 

En reference avec les differents types d'ondes indiquts 
precedemment , la grandeur vectorielle V(M) k estimer peut 
25 etre : 

- un champ eiectrique, dans le cadre de la propagation 
d'une onde eiectrique ; 

- un champ magnet ique, dans le cadre de la propagation 
d'une onde magnetique ; et 

30 - une Vitesse du son au point dans le cadre de la 

propagation d'ondes ultrasonores . 



15 



20 
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Pour determiner les valeurs associSes a chaque source SAi, 
SBi, SCi, on met forme le systeme matriciel selon la 

relation suivante : 



^Vx(Pi) ^ 






Vx(P2) 

• 




VA2 


Vx(Pn) 




vAn 


Vy(Pl) 




vBi 


Vy(P2) 




VB2 


Vy(PN) 




vBn- 


Vz(Pi) 




vCi 


Vz(P2) 

• 
« 




VC2 


.Vz(Pn), 







[11] 



On remargue, en particulier, que la matrice d' interaction 
est de dimensions 3Nx3N, o-Q N est le norabre total 

10 d'Schantillons de surface. La matrice d' interaction 
s'exprime ici par la relation : 

f \ 



C^(i.j) 

cy(i,j) cy(i.j) 

N{3i(i.j) C|(i,j) C^(i,j) 



[12] 
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Les coefficients de cette matrice s'expriment par : 

C:iiiJ)^fu[d(Pi.Sc7j)] tl3] 

5 avec a - A, B, C 

i = 1, 2, N 
j = 1, 2, N 
u = X, y, 2. 

10 et sont encore fonction d'une distance separant le point 
de contact Pi de I'une des sources Soj (o- = A, B ou C) 
d'xin triplet associg a un ichantillon de surface dSj. 

Par inversion de la matrice d' interaction , on 

15 dStermine ainsi les valeurs v<Tj associ§es & chaque source 
So-j, en appliquant des conditions atix limites sur les 
valeurs du vecteur V aux points Pi. Ces conditions aux 

limites imposent une valeur du vecteur V , selon ses trois 
coordonnSes Vx(Pi) , Vy(Pi) et VzCPi) . 

20 

Une fois que ces valeurs de sources vo-j sont ainsi 
dgtermin€es, on calcule facilement 1' expression du vecteur 

V en un point quelconque M de I'espace, par la relation : 
25 V(M) = Vx(M) X + Vy(M) y + V^CM) z [14] 
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<T-A,B,C 



a=A,BjC 



X, y et 2 correspondent a des vecteurs unitaires port6s 
par les axes x, y et z de I'espace tridimensionnel . 

Ainsi, la matrice d' interaction Fy , lorsqu'elle est 
appliqu6e ^ un point quelconque M de I'espace, ne comporte 
finalement que trois lignes associ^es chacune S une 
coordonnee de I'espace x, y ou z. 



Pour diff brents types d'ondes, les valeurs des sources voj, 
sont, comme prgc^demment , \me charge 61ectrique pour une 
onde Slectrique, un flux magnet ique pour une onde 
magnStigue, une vitesse du son pour une onde ultrasonore. 

Plus pr6cis6ment, les coefficients de la matrice 
d" interaction se dSterminent h partir des relations 

[3] , [4] et [5] prec^dentes en pr^cisant toutefois que : 



V(M) = -grad[V(M)] tl5] 

V(M) €tant la grandeur scalaire calculge pr6c6demment par 
25 1* Equation [8] . 
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Ainsi, pour 1 ' estimation d'une grandeur vectorielle V au 
point M et pour les types d'onde cit€s ci-avant a titre 
d'exetnple (electrique, magn§tique et ultrasonore) , les 
coefficients de la matrice d' interaction sont 

5 inversement proportionnels au carr§ d'une distance 
s6parant chague source du point M, alors que pour 
1 ' estimation d'une grandeur scalaire V en un point M de 
I'espace, les coefficients de la matrice d' interaction F 
sont simplement inversement proportionnels & cette 

10 distance. Chaque distance implique I'une des sources d'un 
triplet d'un echantillon de surface et un point M de 
I'espace. La matrice d' interaction Fy comporte alors 3N 

colonnes lorsqu'il s'agit de prendre trois sources par 
Schantillon de surface, alors que la matrice d' interaction 
15 F pour 1» estimation de la grandeur scalaire ne comportait 
que N colonnes puisque seule une source par Schantillon de 
surface etait n^cessaire. - 

Plus gSnSralement , on attribue une source par echantillon 
20 lorsque I'on connait des conditions aux limites pour une 
grandeur scalaire et trois sources par Echantillon lorsque 
I'on connalt plutdt des conditions aux limites pour une 
grandeur vectorielle. 

25 On se refdre maintenant a la figure 3 A pour decrire, ^ 
titre illustratif, une application du procedE selon 
1' invention a 1' estimation d'un potent iel electrique en un 
point M de I'espace tridimensioxmel, situE entre deux 
armatures d'xrn condensateur . Les armatures de ce 

30 condensateur sont portSes ^ des potentiels respectifs VI 
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et V2 . La mise en oeuvre de I'^tape a) consiste d'abord a 
tnailler les surfaces respectives des deux armatures* Dans 
I'exemple represents 3ur la figuxe 5A, on n'a repf§sent6 
que detix mailles pour chacjue armature, simplement ll titre 
5 illustratif. 



10 



L' application de I'^tape b) consiste a mettre en forme le 
systeme matriciel impliquant la matrice d' interaction F et 
le vecteur colonne comportant les valeurs des sources Si a 
S4. La multiplication de ces deux matrices permet 
d'obtenir un vecteur colonne comportant les valeurs du 
potentiel en un ou plusieurs points M de I'espace. 



La mise en ceuvre de I'Stape c) du proc§d6 selon 
15 1' invention consiste & appliquer le systeme matriciel a\ix 
points de contact des hemispheres Pi k P4/ de chaque 



echantillon de surface 
relation suivante : 



dSi k dS4. II en rSsulte la 



20 



^V(Pi)^ 
V(P2) 
V(P3) 

lv(P4)j 



1 



1 



1 



1 



PlSi P1S2 P1S3 P1S4 
1111 



P2S1 P2S3 P2S4 

1 1 1 1 

P3S1 P3S21 P3S3 P3S4 
1111 



P4S1 P4S2 P4S3 P4S4 ^ 



V2 
V3 

Ky4j 



[16] 



avec v(Pi) = vCPa) = Vi 

V(P3) « V(P4) = V2 

Vi = qi, V2 = <l2, V3 = qa/ V4 = q4 



25 
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Ici, la condition aux limites impose que la valeur du 
potentiel aux points de contact Pi et P2 correspond au 
potentiel VI de la preniSre armature, De meme, le 
potentiel electrique axix points de contact P3 et P4 
5 correspond au potentiel Electrique de la seconde arinature 
V2, Par inversion de la matrice d ' interaction appliqu^e au 
point de contact Pi, on determine les valeurs des sources 
Vi qui correspondent, comme exprimg dans la relation [16] , 
a des charges Electriques q±. 

10 

Les coefficients de la matrice d' interaction sont 

parfaitement connus, puisque les positions des sources Sj 
et les positions des points de contact Pi sont dSterminSes 
au pr^alable, comme represents sur la figure IB. 

15 

L' expression du potentiel Electrique V(M) au point M entre 
les deux aanmatures est finalement donnge par 
1 ' es^ression : 



20 V(M) = ^fjil- + ^ + ^+ ^4 ^ 



Ittsq \^MSi MS2 MS3 MS4 J 



[17] 



On se rEfSre maintenant §l la figure 3B sur laquelle on a 
represent E les m§mes armatures que sur la figure 3A, avec 
sensiblement un mSme maillage, mais dans le but, ici, 
25 d'estimer \me grandeur vectorielle correspondant au cliamp 

Electrique E(M) , au point M de I'espace tridimensionnel . 
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On peut appliquer les relations [11] k [15] pour estimer 
la valeur du champ electrique au point M, avec, dans la 
relation [3.3] : 



[18] 



avec 



10 



cr = A, B, C 

i = 1, 2, 3, 4 
j = 1, 2, 3, 4. 



Toutefois, les valeurs du cliamp Electrique au point de 
contact Pi restent k determiner dans la relation [11] . 

15 On introduit alors une loi ggnerale pr6d6termin§e du 
comportement du cliamp (en reflexion, en transmission, ou 
autre) au niveau de la surface de 1" obstacle (des 
armatures dans I'exemple du condensateur precitE) , pour 
connaitre les valeurs des sources Vaj . 

20 

Par exemple, si l»onde electrique est totalement reflechie 
par la surface d»un obstacle (par exemple I'une des deux 
armatures) , le cliamp Electrique en un point de contact Pi 
est normal la surface dSi et ses composantes Ex et Ey 
25 sont nulles. A titre illustratif, si la surface de 
1» armature n'Etait reprSsentSe que par un seul Echantlllon 
de surface k trois sources, les valeurs de ses sources vA, 
vB et vC seraient toutes Egales entre elles k une mime 
valeur +q. 
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Au contraire, si le coefficient de reflexion est 
pratiquement nul k la surface dSi, la composante du champ 
Slectrique au point Pi est nulle, ce qui correspond 
5 bien au cas oil le champ est sensiblement tangent H la 
surface dSi. Ainsi, a titre illustratif, si la surface de 
1' armature n'gtait representee que par un seul gchantillon 
de surface a trois sources, les valeurs de ses sources vA, 
vB et vC seraient respectivement , par exemple, +q, +q et 

lb -2q. Par exemple, dans le cadre de la propagation d'une 
onde raagnetique, si la surface d'un capteur a courants de 
Foucault (avec une composante normale du champ magn§tique 
nulle) etait representee par un unique echantillon de 
surface, les flux magnetiques des trois sources associees 

15 a cet echantillon de surface seraient +(p,+q>et -2<p. 

On comprend ainsi qu'avec les trois sources par 
echantillon dSi, on peut definir, par exemple en fonction 
de la ponderation de chaque source, une orientation 
20 quelconque du champ t la surface de 1' obstacle. 

Bien entendu, cette demarche suppose que le coefficient de 
reflexion R d'xin obstacle soit connu auparavant. En 
particulier, il peut itre avantageiox de comparer une 

25 simulation et une mesure experimentale pour detecter, a la 
surface d'un obstacle, des inhomogeneites ou des irapuretes 
qui correspondent a des points de la surface de cet 
obstacle qui ne verifient pas les valeurs du coefficient 
de reflexion R imposees a chaque point predetermine Pi de 

30 1 'obstacle. 
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On peut affecter ainsi \me valeur prSdSterminSe du 
coefficient de reflexion k chaque point Pi de la surface 
de 1' obstacle. A cet effet, oxx incroduit line matrice R qui 
est representative du coefficient de reflexion en chaque 

5 point Pi. Pour une interaction entre un element rayonnant 
et un obstacle, on peut ainsi exprimer le systeme 
matriciel de la relation [9] autrement, c'est-a-dire en 
donnant une expression unique de toutes les sources du 
systSme {k la fois de 1' obstacle et de 1' element 

10 rayonnant), coirane indiqu6 ci-aprSs. 



Dans ce qui suit, on indique que : 

- P(P) est la matrice d' interaction de 1' obstacle OBS 
appliqu€e aux points Pi de la surface de 1' obstacle 

15 OBS ; 

- P(P') est la matrice d' interaction de 1' obstacle OBS 
appliquee aux points P'l de la surface de 1» element 
rayonnant ER; 

FMP) correspond a la matrice d' interaction de 
20 I'^l^ment rayonnant ER appliquee aux points Pi de la 

surface de 1 ' obstacle OBS ; 

- F"(P») correspond & la matrice d» interaction de 
1' Element rayonnant ER appliquee aux points Pi de la 
surface de 1 ' Pigment rayonnant ER ; 

25 - v' correspond au vecteur colonne comportant les valeurs 
des sources S ■ i de 1 ' 61€ment rayonnant ER ; et 
V correspond au vecteur colonne comportant les valeurs 
des sources Si de 1' obstacle OBS. 
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Sur 1' obstacle, la contribution de I'onde gmise par 
I'SlSment rayonnant ER s'exprime par : 

V'(P) = P(P).v' [193 

5 

La contribution de I'onde secondaire renvoy^e par 

1» obstacle OBS s'exprime, par definition, par la 
relation : 

10 V(P) = F(P).v [20] 

Or, dans 1 ' example repr^sente sur la figure 4A, I'onde 
secondaire correspond simplement a une reflexion de I'onde 
principale, Ce qui s'exprime par la relation : 



15 



V(P) = RV'(P) [21] 



oil R correspond une matrice de reflexion dont chaque 
coefficient represente la contribution k 1' Emission, par 
20 reflexion, de I'onde secondaire, par chaque source Si (ou 
SOi, dans le cadre d'une estimation d'une grandeur 
vectorielle) de 1' obstacle OBS. 



Des trois relations [19] , [20] et [21] , on d^duit 
25 1' expression du vecteur colonne v comportant les valeurs 
des sources de 1' obstacle, t partir du vecteur colonne v 
comportant les valeurs des sources de l'61#ment rayonnant, 
par la relation : 



30 v^[FiP)]'^ .R.[F'(P)].r 



[22] 
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Par ailleurs, pour \ine estimation fine des grandeurs 
scalaires ou vectorielles au poini. hi, nocammeni: pour tenir 
compte de reflexions multiples, il est avantageux de tenir 

5 compte de la contribution du rayonnement par 1' obstacle, 
au niveau de la surface de I'^l^ment rayonnant ER. A cet 
effet, on tient compte, dans 1' estimation des conditions 
aoax limites a la surface de 1' element rayonnant ER (aux 
points P'i) de la contribution du rayonnement des sources 

10 S'i de 1' element rayonnant et de la contribution de 
1» Emission de I'onde secondaire par les sources Si de 
1 'obstacle, par la relation : 

Vt(F) = F(F) V + P(F) v* [23] 

15 

On peut ainsi ajuster, grSce a la relation [23] , les 
valeurs de sources S'l de 1' Element rayonnant ER, en 
tenant compte de la reflexion de 1 ' obstacle OBS, selon la 
relation suivante : 

20 

Vt(P') = ^(P') .[F(P) .R .[F'(P) ] + F'(P') }.v' [10] 

Ainsi, on impose simplement des conditions aux limites 
pour 1' Element rayonnant, pour en d€duire les valeurs des 
25 sources v'l. En pratique, on procedera pref €rentiellement 
comme suit : 

- apres maillage des surfaces, on determine la position 
des points Pi et P'l et des sources Si et S'l; 

- en fonction du type d'onde en jeu, on determine les 
30 coefficients des matrices P(P) , F' (P) , F(P') et F'(P'); 
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- en fonction d'une loi de rgflexion de 1 'obstacle, on 
dgtermine les coefficients de la matrice de r6flexion 
comme dans i'epperople donnS plus loin pour une onde 
ultrasonore ; 

- en fonction de conditions aiax limites sur 1' Element 
rayonnant (dont le comportement est ggn^ralement connu 
pour un probldme donn€) , on determine les valeurs du 
yecteur ViiP^ aux points P'i de la surface de 1- Pigment 
rayonnant et on en diduit les valeurs des sources S'i de 
1" Pigment rayonnant par inversion de la relation [10]; 

- on en dSduit aussi les valeurs des sources Si de 
1' obstacle par application de la relation [22]; 

- une fois les valeurs de toutes les sources Si et S'i 
dSterminSes, on peut appliquer le systdme matriciel donng 

15 par la relation [9] a tout point M de I'espace, en 
appliquant k ce point M les matrices d« interaction F et P' 
(impliquant la position du point M et les positions des 
sources respectives Si et S'i). 



10 



20 



25 



En se r4fgrant a nouveau a la figure 4A, on considdre que 
1' obstacle OBS reprgsente simplement une interface entre 
deux milieux Ml et M2, formant ainsi un dioptre qui peut 
gtre plan, tel que repr^sente dans I'exemple de la figure 
4A, mais aussi incurvg ou de forme ggn€rale quelconque. 
Les coefficients de reflexion Ri associgs k chaque point 
Pi dependent, dans le cadre de la propagation d'uxie onde 
ultrasonore ou €lectromagn#tique de haute frequence, de 
1' angle d' incidence pi du rayon issu de la source Si, au 
point de I'espace tridimensionnel M. 
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Pour une onde ultrasonore, 1' expression des coefficients 
de reflexion Ri est donn6e par : 



Ri = 



/?2C2 COSfii - piCi 


2 2 

+-2" COS Pi 

. cf of 


1/2 


P2C2 COSyffi + PiCi 


l-4^4eos2^/ 
- cf cf 


1/2 



[24] 



10 



15 



10 



- ci est la Vitesse du son dans le milieu Mi ; 

- C2 est la Vitesse du son dans le milieu M2 ; 

- pi est la density du milieu Mi ; 

- P2 est la density du milieu M2. 

Dans cette expression [24] ; le terme cosPi peut simplement 
§tre estimg en fonction des coordonn^es dans I'espace du 
point M et du point reprgsentant la source Si. 

En se r#f grant maintenant a la figure 4B, la mgme 
estimation peut gtre menee pour xin point M situS dans le 
milieu M2. Dans ce cas, I'onde que re<?oit le point M est 
une onde transmise par 1' obstacle OBS. En particulier, on 
remarque que les sources de 1' Pigment rayonnant ER ne sont 
plus actives, du fait de 1 'occultation de 1' Pigment 
rayonnant ER par 1- obstacle OBS. En transmission, le 
raisonnement s' applique comme pr6c6demment avec \ane 
condition aux limites impos6e aux points Pi par les 
valeurs des coefficients de transmission Ti associ#s a 
chaque point Pi. Dans le cadre de la propagation d'une 
onde ultrasonore, chaque coefficient de transmission Ti 
est donnS par la relation : 
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Ti = 



1/2 



[25] 



of 



5 Comme indiquS prScSdemment , les termes cospi peuvent @tre 
d#termin6s en fonction des coordonnees respectives des 
sources Si et du point M. 

Pour estimer les valeurs de sources Si de 1' obstacle OBS, 
10 on applique la relation [22] en remplagant toutefois la 
matrice de reflexion R par la matrice de transmission T : 



v-[F(P)]"^T[P(P)]v' [26] 

15 Dans le cadre de la propagation d'une onde ultrasonore, 
les coefficients des matrices R et T s'estiment pour 
ctiaque source Si et pour cliaque point Pi. En particulier, 
cliaque coefficient Ti,j ou Ri,j de la matrice T ou de la 
matrice R (od i correspond k la i^""® ligne et j correspond 

20 a la j^"^® colonne) s'exprime en fonction d"un angle fi±^ 
entre line normale ^ la surface de 1 ' obstacle au point Pi 
et une droite passant par le point Pi et par line source 
Sj. On peut Scrire ainsi, de manidre gSnSrale, les deux 
relations exprimant les valeurs des coefficients des 

25 matrices T et R par les relations respectives suivantes : 



= ft(C0S)5ij) 



[27] 
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Ry = fr(cos^y) [28] 

Oil ft est donn6 par la relation [25] et fr est donn§ par 
la relation [24] . 

5 

De fagon plus g€n6rale, en r6f€rence aiix figures 4A et 4B, 
on indique que, si I'on considere 1* obstacle conunae un 
mat^riau plein reprSsentant un milieu M2 distinct d'un 
milieu Ml dans lequel se propageait initialement I'onde 

10 principale : 

- pour tine reflexion de I'onde principale sur 1" obstacle 
en tant que milieu M2 (la surface de 1' obstacle formant 
dioptre entre les milieux Ml et M2) , les ligmisptidres HEMi 
sont orientes vers 1 ' exterieur de 1' obstacle (figure 4A) ; 

15 - pour une transmission de I'onde principale dans 
1' obstacle, les li^mxspheres HEMi sont orientes vers 
l»int6rieur de 1' obstacle (figure 4B) . 



On se r€fdre maintenant k la figure 5A pour d6crire le cas 
20 d'un 1 'obstacle plan OBS de dimensions finies, excite par 
un element rayonnant ER, incline d'un angle pr6d6termin6 
par rapport t 1' obstacle OBS. Comme indiquS pr6c§demment , 
pour une onde ultrasonore, on prendra en compte 
1' inclinaison de I'SlSment rayonnant pour calculer la 
25 contribution de I'onde 6mise par l'61§ment rayonnant au 
point M. Par ailleurs, de fagon particuliSrement 
avantageuse, on maille une surface qui englobe la surface 
de 1» obstacle (figure 5A) . Pour une tranche de I'espace 
delimit^e par 1' element rayonnant, d'une part, et 
30 1» obstacle, d' autre part (figure 5A) , on peut considerer 
trois types de sources : 
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les sources S*± de 1» element rayonnant ER, 
des sources SOi, qui renvoient I'onde secondaire, par 
reflexion de 1! obstacle OBS, en fonction d'un certain 
coefficient de reflexion R de 1' obstacle ; et 
des sources S3±, qui ne renvoient pas d'onde secondaire 
et auxquelles un coefficient de reflexion nul peut etre 
attribug si 1' obstacle sepaire deux milieux d' indices 
identiques. Dans ce cas, ces sources SSi sont 
consid€r§es comma "Steintes" dans la tranche de 
I'espace prScitSe et ne sont pas prises en compte dans 
les calculs de la grandeur physique au point M de la 
figure 5A. En revanche, ces sources SSi peuvent §tre 
actives par reflexion de I'onde principale si 
1' obstacle OBS s6pare deux milieux d' indices 
dif f grents . 

Par ailleurs, pour estimer les grandeurs scalaires ou 
vectorielles associees a un point M d'un demi-espace 
delimite par la surface englobant 1' obstacle OBS (a droite 
. 20 de la figure 5B) , on considdre : 

- les sources SO'i de 1' obstacle, actives par transmission 
de I'onde principale, et 

les sources SSi, auxquelles on affecte maintenant un 
coefficient de transmission Sgal ^ 1 si 1* obstacle 
25 sSpare deux milievix de m€mes indices. Ces sources SSi se 

comportent finalement (aux angles d' incidence prds) 
commg les sources S*i de 1' element rayonnant ER. 
Les sources S'l de 1' Pigment rayonnant peuvent alors §tre 
"6teintes" pour le calcul des grandeurs physiques dans ce 
30 demi - e space . 



10 
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Pour calculer les valeurs v* des sources S'i de I'Sl^ment 
rayonnant (desquelles sont dSduites les valeurs v des 
sources de 1' obstacle selon les xelafcions [22] et [26J ) , 
on appliquera s implement des conditions aux limites aux 
5 points de la surface active de 1' element rayonnant ER. Par 
exemple, pour une propagation d'onde ultrasonore, on peut 
indiquer cjue les vitesses sonores aux points de la surface 
de 1' element rayonnant ER sont perpendiculaires k cette 
surface et de modules vo egaux entre eux. 

10 

De maniSre g#n#rale, on indique que I'espace 
tridimensionnel peut §tre ainsi decoupe par des interfaces 
dSlimitant des milieux de proprietSs distinctes, chaque 
interface reprgsentant un obstacle au sens de la prSsente 

15 invention. On calcule alors les grandeurs physiques dans 
chaque tranche de I'espace. Par exemple, dans le cadre de 
l»6tude d'une hStSrostructure {k plusieurs interfaces), on 
peut appliquer le precede ci-avant pour des tranches 
successives de 1 ' espace en considSrant deux interfaces : 

20 I'une reprSsentant un "Element rayonnant" au sens des 
figures 4A et 5A, par exemple par transmission d'une onde 
regue, et 1' autre repr^sentant un obstacle recevant I'onde 
transmise. Avantageusement , on tient compte, pour chaque 
tranche de I'espace, des contributions de toutes les 

25 interfaces comme exprime par les relations [10] et [22] . 

Toutefois, dans xine realisation pratique prSf^rSe, en 
particulier pour programmer la simulation d'une 
interaction, on considgrera avantageusement tous les 
30 obstacles de tout I'espace autour d'un point M et I'on 
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itnposera une condition sur la position du point M par 
rapport k chaque source prSsente dans I'espace. 

Pr^f^rentiellement, en rif^rence & la figure 7A, on teste 

5 le produit scalaire SMj- k chaque iteration par rapport a 
une source S, par exenqple sous la forme : 

rii]Hr'°' 

oil r est le vecteur reliant la source S au point de 
contact P de la demie sphere avec I'^l^ment de surface dS 
10 considers, dans le cas oil I'on privoit une seule source 
par hSmisphSre. Dans le cas ovl I'on prgvoit plut6t un 

triplet de sources SI, S2, S3, la base du vecteur r est 
pr^f^rentiellement situee au barycentre des trois sources 
SI, S2, S3. D'ailleurs, dans le cas de trois sources par 
15 gchantillon de surface, le calcul du produit scalaire 
conceme chaque source Si du triplet SI, S2, S3. 

Typiquement, si cos^ est positif, on tient compte d^ la 
contribution de cette source S dans 1 ' estimation de 
20 1 ' interaction . 

En revanche, si cos^ est n^gatif, on affecte une valeur 
(scalaire ou vectorielle) nulle k cette source S dans 
1 ' estimation de 1 • interaction. 

25 

En variante du calcul du produit scalaire ci-avant, on 
peut calculer une "altitude" du point M. Dans ce cas, le 
test porte sur une grandeur du type : 
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SM.r 

ihi 

Bien entendu, d'autres types de tests sont possibles. Par 
5 exemple, dans le cas d'lm calcul par rapport a 1' angle 0, 
on peut choisir cet angle dans un c6ne d"ouverture 
choisie, ou autre. 

Pinalement, cette d§marche pertnet avantageusement de 
10 systSmatiser ime configuration quelconque des sources par 
rapport au point d' observation M, en introduisant 
simplement une gtape supplement a ire de test, k chaque 
iteration sur une source de la position de cette source 
S par rapport au point M, comme indiqu§ ci-avant. 

15 

Cette demarche s'avere particulierement avantageuse pour 
des surfaces mailler qui sont relativement complexes, 
notamment lorsque le point d' observation M est susceptible 
de se situer dans \me zone d' ombre par rapport k certaines 

20 sources, comme le montre la figure 7B. Sur cette figure 
7B, on a repr§sent# en traits pointings les demi- spheres 
associSes aux sources dans la zone d' ombre du point 
d' observation M, et pour lesquelles^ par consequent, la 
contribution est . fix6e comme 6tant nulle dans 

25 1 » interaction estimee. Bien entendu, dans le cas oOl une 
source est 6crantSe par une surface 6chantillonn6e, bien 

que le produit scalaire SMj* associe k cette source reste 

positif, un second test determine si le vecteur SM 
traverse ou non vine surface Schantillonn^e. Dans 
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1' affirmative, cette source est consid#ree comme inactive 
specif iquement pour la rSgion du point M. 

Ainsi, en termes plus ggn^raux, le precede au sens de 
5 1 » invention prevoit pref erentiellement au moins une etape 
supplement aire, pour chaque §chantillon de surface, de 
test de la valeur d'un produit scalaire entre : 

- un premier vecteur r normal k 1 ' 6chantillon de surface 
et dirigS vers le sommet P de 1 'hSmisphSre,. tel que 

10 represents sur la figure 7A, et 

- un second vecteur SM tir# entre une source S associSe a 
cet hemisphere et le point M qui se situe dans la region 
d' observation, 

en distinguant, en particulier : 
15 - le cas oil ce produit scalaire est infSrieur a un 

seuil prSdSterminS et la contribution de cette source 
n'est pas prise en compte, et 
- le cas ou ce produit scalaire est supSrieur a un 
seuil predetermine et la contribution de cette source 
20 est ef f ectivement prise en compte. 

Dans I'exemple ci-avant oil I'on considere 1» angle 0 entre 
ces deux vecteurs, le seuil predetermine precite est bien 
entendu la valeur 0 et I'on distingue simplement les cas 
25 oia le produit scalaire est positif ou negatif • 

Bien entendu, ce choix n^est pas limitatif de sorte que 
pour line heterostructure k plusieurs dioptres paralieies, 
on pourra encore considerer, de fagon avantageuse, des 
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semi-espaces successifs, cotnme d^crit ci-avant en 
r§f§rence avuc figures 4A et 4B. 

La simulation de la figure 5C correspond, pour une onde 
ultrasonore, h la situation des figures 5A et 5B en tenant 
corapte : 

- de la contribution de 1- Emission de I'onde principale 
par I'SlSment rayonnant ER ; 

- de la contribution de la reflexion de cette onde 
principale par 1« obstacle ; et 

- de la contribution de la transmission de I'onde 
principale par 1' obstacle. 

Les lignes de niveau de la figure 5C correspondent a 
diff€rents paliers de pression acoustique. L'616ment 
rayonnant ER est dispose a 10 mm de 1' obstacle OBS et 
incline de 20^ par rapport a ce dernier. On remarque en 
particulier des f ranges d- interferences dans une zone 
entre 1' obstacle OBS et l'61€ment rayonnant ER. Une telle 
simulation peut avantageusement indiquer une position 
id^ale d'un capteur ultrasonore. Ces capteurs ultrasonores 
comportent habituellement un transducteur comme 616ment 
rayonnant actif et un d^tecteur pour mesurer les ondes 
ultrasonores regues. La simulation de la figure 5C peut 
ainsi indiquer en outre la forme id^ale d'un capteur 
ultrasonore, selon les applications souhait6es, pour une 
forme d' obstacle donnSe. 

La simulation de la figure 5C a M ef fectu^e grace un 
calcul matriciel programme a I'aide du logiciel de calcul 
MATLAB©. Le nombre de mailles total choisi pour 
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1' obstacle et pour 1' element rayonnant (ici^ quelques 
centaines en tout) est alors optimise : 

- d'une part, pour limiter la durfee des calculs ; et 

- d' autre part, pour que la taille des mailles reste 
5 inf^rieure & une demie longueur d'onde, de maniere S 

verifier le crit^re de Rayleigh. 

On indique toutefois que, corame les ^ISments & mailler 
dans la mise en oeuvre du proc€d6 selon 1' invention sont 
10 simplement des surfaces, les temps de calcul sont bien 
moins longs que ceux n€cessaires dans la tnise en oeuvre 
d'un proc#d§ de calcul de type "par ^I^ents finis". 

La prSsente invention peut ainsi se manif ester par la mise 
15 en oeuvre d'xane succession d' instructions d'un produit 
programme d'ordinateur stocke dans la mgmoire d*un disque 
dur ou sur un support amovible et se d^roulant comme 
suit : 

choix d*un pas de maillage notamment en fonction de la 
20 longueur d*onde de I'onde principale ; 

- determination des coordonn^es des sources Si et/ou S ' ± 
et des points de contact Pi et/ou P'i ; 

- choix d'un type d'onde en jeu et calcul des 
coefficients des matrices d' interaction appliqu6es axix 

25 points Pi et/ou P'i par la mise en oeuvre d'un logiciel 

de calcul matriciel ; 

- choix d'une loi de reflexion et/ou de transmission de 
la surface de 1' obstacle et calcul des coefficients des 
matrices de reflexion et/ou de transmission ; 

30 - calcul des valeurs des sources Si et/ou S'i ; et 
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- calcul des grandeurs scalaires ou vectorielles en tout 
point de I'espace tridimensionnel . 

A ce titre, la pr^sente invention vise aussi un tel 
produit programme d'ordinateur, stocks dans une mSmoire 
d' unite centrale ou sur un support amovible propre a 
coop§rer avec tin lecteur de cette unit§ centrale, et 
coraportant en particulier des instructions pour mettre en 
OBUvre le proc6d6 selon 1' invention. 

Bien entendu, la pr^sente invention ne se limite pas k la 
forme de realisation dScrite ci-avant & titre d'exeraple ; 
elle s'etend k d'autres variantes. 

Ainsi, on comprendra que, m§me si dans les figures 
cotnmentees ci-avant on repr^sente a la fois la surface 
d'un obstacle et la surface d'un element rayonnant, la 
pr^sente invention s' applique aussi & 1' estimation de 
grandeurs physiques dans le cadre d'une onde interagissant 
avec vin obstacle et 6mise en champ lointain. Dans ce 
contexte, il n«est pas n^cessaire de matSrialiser la 
surface d'un Aliment rayonnant H mailler et les relations 

[8] et [14] ci-avant suffisent a determiner 1 • interaction 

entre cette onde et 1« obstacle. 

On a indique ci-avant des Equations perroettant de calculer 
les grandeurs scalaires ou vectorielles en un point M de 
I'espace, pour des ondes glectromagngtiques, ou encore 
acoustiques. Bien entendu, ces grandeurs peuvent Stre 
estimges pour d'autres types d' ondes, notamment pour des 
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ondes therraiques, des ondes 61ectromagn6tiques impliquant 
des antennes radiof r^tjuences , ou autres. 

Bien entendu, la pr^sente invention ne se limite pas k xine 
5 application au contrdle non destructif, mais a tout type 
d' application, notatnment en imagerie m^dicale, par example 
pour l'6tude de microsystdmes mettant en ceuvre une 
microscopie acoustique & miroirs mobiles. 

10 On a dScrit ci-avant des interactions entre \ine onde et un 
obstacle unique. Bien entendu, la pr^sente invention 
s' applique a une interaction avec plusieurs obstacles. A 
cet effet, il convient simplement de mailler les surfaces 
de ces obstacles et d' additionner leur contribution pour 

15 1' estimation d'\ane grandeur vectorielle ou scalaire en \in 
point quelconque de I'espace. De m§me, comme indiquS ci- 
avant, la surface de 1' obstacle OBS peut Stre plane, ou 
encore incurv^e, ou encore de forme complexe quelconque. 

20 Ainsi, dans le cadre d'une onde interagissant avec 
plusieurs obstacles dans I'espace, \ine simulation 
gquivalente a celle representee sur la figure 5C 
permettrait de positionner des capteurs et/ou des elements 
rayonnants en fonction de la configuration de ces 

25 obstacles, notamment pour \ine application a la 
determination de la position de hauts-parleurs dans un 
habitacle cloisonne, tel qu'un habitacle de vehicule 
automobile . 

30 L'espace tridimensionnel peut etre decoupe en une 
pluralite de regions, comme decrit ci-avant en reference 
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atax figures 4A, 4B, 5A et 5B. Toutefois, pour qu'\ane 
surface de I'une desdites regions soit considSrSe comme un 
obstacle d'une onde prinuipale, devenant actif par 
emission d • une onde secondaire , 1 ' incidence de 1 • onde 
5 principale sur cette surface doit pr6f ^rentiellement 
rester inf§rieure ou ggale & 90°. 
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Revendications 



1. Proc6d6 ppur evaluer une grandeur physique associSe ^ 
une interaction entre une onde et un obstacle, dans une 
5 region de 1 • espace tridimensionnel, dans lequel ; 

a) on determine, par maillage, une plurality 
d'echanti lions de surface (dSi) dont une partie au moins 
repr^sente la surface d'un obstacle recevant \xne onde 
principale et §mettant, en r§ponse, xme onde secondaire, 

10 et I'on attribue ^ chaque §chantillon de surface au moins 
une source (Si) Smettant une onde Element aire reprSsentant 
une contribution k ladite onde secondaire, 

b) on forme un systSme matriciel coraportant : 

- une matrice d' interaction {P(M)), inversible, 
15 appliquSe & une region donnSe (M) de 1' espace et 

comportant un nombre de colonnes correspondant ^ un 
nombre total de sources, 

une premiere matrice colonne dont chaque coefficient 
(vi) est associe a une source (Si) et caracterise 

20 I'onde glementaire qu'elle 6met, 

et une seconde matrice colonne, obtenue par une 
multiplication de la premiSre matrice colonne par la 
matrice d' interaction et dont les coefficients sont 
des valeurs d*une grandeur physique (V(M)) 

25 representative de I'onde emise par 1' ensemble des 

sources en ladite region donnee (M) , 

c) pour estimer les coefficients de la premiSre matrice 
colonne (vi) , on af facte des valeurs de grandeur physicjue 
choisies {V(Pi)) k des points predetermines (Pi), associ^s 

30 chacun k un ^chantillon de surface (dSi) , lesdites valeurs 
choisies (V(Pi)) 6tant rangSes dans la seconde matrice 
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colonne, et I'on multiplie cette seconde matrice colonne 
par 1' inverse de la matrice d' interaction appliqu^e 
auxdits points pr^d^terminSs (Fi) , 

d) pour ^valuer ladite grandeur physique (V(M)) 
5 reprSsentant I'onde Smise par 1' ensemble des sources dans 
une region donn6e (M) de I'espace tridimensionnel , on 
applique la matrice d' interaction k ladite region donn€e 
(M) et on multiplie cette matrice d' interaction par la 
premiere matrice colonne comportant les coefficients 
10 estim6s & l'6tape c) . 

2. Procgdg selon la revendication 1, dans lequel, pour 
^valuer une grandeur physique representative d'une 
interaction entre un Pigment rayonnant une onde principale 
15 et un obstacle recevant cette onde principale, 

- a I'gtape a), on determine en outre, par maillage, une 
plurality d>€chantillons de surface (dS'i) reprSsentant 
ensemble une surface active de I'^l^ment rayonnant I'onde 
principale, et I'on attribue a chaque gchantillon de la 

20 surface active au moins une source (S'i) ^mettant une onde 
glgmentaire reprSsentant une contribution k ladite onde 
principale, 

- on applique en outre les Stapes b) , c) et d) aux 
gchantillons de la surface active, et 

25 - on 6value ladite grandeur physique (V(M)) reprSsentant 
1 ' interaction entre 1' Pigment rayonnant et 1' obstacle dans 
une region donnSe (M) de I'espace tridimensionnel, en 
prenant en compte la contribution, en ladite region donnSe 
(M), de I'onde principale 6mise par 1' ensemble des sources 

30 de la surface active et la contribution de I'onde 
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secondaire Smise par 1' ensemble des sources de la sijirface 
de 1' obstacle. 

3. Proc6d^ selon I'une des revendications 1 et 2, dans 
5 lequel chaque coefficient de la matrice d' interaction, 

appliqu6e & une region donnee de I'espace, est 
repr^sentatif d'une interaction entre une source et ladite 
region donnee et la valeur de chaque coefficient est 
fonction d'une distance entre une source et ladite r€gion 
10 donn6e • 

4. ProcedS selon I'une des revendications 1^3, dans 
lequel la matrice d' interaction appliquSe, k I'etape c) , 
aiixdits points pr€d6termin6s (Pi) , comporte un nombre de 

15 lignes correspondant k un nombre total de points 
predStearminSs (Pi) . 

5. Proc^dg selon I'une des revendications 1 ^ 4, dans 
lequel la grandeur physique a ^valuer est une grandeur 

20 scalaire (V(Pi)) et, a I'etape a), on attribue oine source 
unique k chaque echantillon de surface. 

6. ProcSdS selon la revendication 5, dans lequel la 
matrice d" interaction (F (M) ) appliqu^e, k 1 ' Stape d) , k 

25 une region de I'espace (M) , comporte \me ligne. 

?• Proc€d6 selon I'une des revendications 5 et 6, dans 
lequel chaque point pr#d6termin6 (Pi) associS k un 
Echantillon de surface (dSi) correspond k un point de 
30 contact entre cet Echantillon de surface (dSi) et un 
hEmisphSre dont la surface est Egale k la surface de cet 
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gchantillon de surface, et de centre correspondant a une 
position de la source (Si) qui est attribute k cet 
Schantillon de surface. 

5 8. Proc6d6 salon I'une des revendications 5 k 7, dans 
lequel : 

- I'onde principale est une onde §lectrique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de charge glectrique associ^es chacune k \ine 

10 source, et 

- les coefficients de la seconde matrice colonne sont des 
valeurs de potentiel ^lectrique, 

9. ProcSdS selon I'une des revendications 5^7, dans 
15 lequel : 

- I'onde principale est une onde magnStique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de flux magnStique associ^es chacune Sl une source, 
et 

20 - les coefficients de la seconde matrice colonne sont des 
valeurs de potentiel magnetique. 

10. Precede selon I'une des revendications 5^7, dans 
lequel : 

25 - I'onde principale est une onde sonore, 

- les coefficients de la premidre matrice colonne sont des 
valeurs de vitesse de son associ^es chacune k une source, 
et 

- les coefficients de la seconde matrice colonne sont des 
30 valeurs de pression acoustique. 
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11. Proc€di selon I'une des revendi cat ions 1^4, dans 
lequel la grandeur physique Bl ^valuer est une grandeur 
vectorielle (V(Pi)) exprimee par ses trois coordonn#es 
dans I'espace tridimensionnel, et I'on attribue, & l'6tape 
a), trois sources (SAi, SBi, SCi) a chaque gchantillon de 
surface (dSi) . 

12. Proc6d6 selon la revendication 11, dans lequel la 
matrice d • interaction (Fv{M)) appliqu^e, k l'§tape d) , a 
\ine region de 1 ' espace (M) , comporte une ligne pour chaque 
coordonn^e (X,Y,Z) de I'espace. 

13. Proc€d6 selon l'\ine des revendi cat ions 11 et 12, dans 
lequel : 

- les trois sources attributes k chaque 6chantillon de 
surface sont sensiblement dans un m^me plan et 

- chaque point predetermine (Pi) associt k un gchantillon 
de surface (dSi) correspond & un point de contact entre 
cet gchantillon et un hemisphere dont la surface est ggale 
a la surface de cet echantillon, et de centre 
correspondant k la position d'un barycentre des trois 
sources . 

14. ProcSdS selon la revendication 13, dans lequel les 
trois sources d'un mSme tchantillon de surface forment 
sensiblement un triangle Equilateral, et les triangles des 
echantillons de surface sont orientes sensiblement 
aieatoirement les uns par rapport aux autres. 

15. Precede selon I'une des revendications 11 k 14, dans 
lequel : 



15 



20 



wo 2004/044790 



PCT/FR2003/003323 



52 

- I'onde principals est une onde Slectrique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de charge Slectrique associ^es chacune & une 
source , et 

- les coefficients de la seconde matrice colonne sont des 
valeurs de cooirdonnSes d'lin champ ^lectrique. 

16. Proc#d6 selon I'vme des revendications 11 a 14, dans 
lequel : 

- l«onde principale est une onde magnet ique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de fl\ix magnet ique associSes chacune a une source, 
et 

- les coefficients de la seconde matrice colonne sont des 
valeurs de coordonn^es d'un chan^ magnStique. 

17. Proc6di selon I'une des revendications 11 k 14, dans 
lequel : 

- I'onde principals est une onde sonore, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de vitesse de son associ^es chacune a une source, 
et 

- les coefficients de la seconde matrice colonne sont des 
valeurs de coordonn^es d'xine vitesse acoustique. 

18. ProcSdS selon I'une des revendications pr6c6dentes, 
dans lequel, pour estimer la contribution de I'onde 
secondaire en ladite region donnie Sl I'Stape d) , lesdites 
valeurs de grandeur physique (V(Pi)) choisies k l'6tape c) 
sont fonction d'un coefficient pr6d§termin6 de reflexion 



10 



15 



wo 2004/044790 



PCT/FR2003/003323 



53 

et/ou de transmission de l^onde principals par chaque 
§chantillon de surface de 1' obstacle. 

19. Precede selon la revendication 18, prise en 
5 cotnbinaison avec I'une des revendications 6 et 12, dans 
lequel I'onde secondaire correspond k une reflexion de 
1 ' onde principale sur 1 ' obstacle et 1 ' hemisphere est 
oriente vers 1 * exterieur de 1 ' obstacle . 

lb 20. Proc6d6 selon la revendication 18, prise en 
corabinaison avec I'une des revendications 6 et 12, dans 
lequel I'onde secondaire correspond a une transmission de 
1 • onde principale dans 1 ' obstacle et 1 ' h€mi sphere est 
orients vers 1 • intSrieur de 1 • obstacle . 

15 

21. Proc6d6 selon I'tine des revendications 19 et 20, dans 
lequel, k I'Stape c) , on determine les valeurs (v'i) 
associSes aux sources (S'l) de 1" Pigment rayonnant (ER) et 
on met en forme au moins : 
20 ~ une premiere matrice d' interaction (F(P)) representant 
la contribution des sources de 1' obstacle aux points 
predetermines de la surface de 1' obstacle (Pi), 

- une seconde matrice d' interaction (F'(P)) representant 
la contribution des sources de 1« element rayonnant aux 

25 points predetermines de la surface de 1» obstacle (Pi), 

- une matrice de reflexion (R) ou de transmission (T) , 
dont les coefficients reprSsentent des coefficients de 
reflexion ou de transmission en chaque point predetermine 
(Pi) de 1' obstacle, 

30 pour determiner les valeurs des sources de 1' obstacle (vi) 
en fonction des valeurs des sources de 1« element rayonnant 
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(v'i) et d'une multiplication des premidre et seconde 
matrices d' interaction et de la matrice de rSflexion ou de 
transmission. 

5 22. Proc^dg selon la revendication 21, dans lequel, a 
I'fitape c) , on determine les valeurs (v'i) associ^es aux 
sources (S'i) de I'^l^ment rayonnant (ER) en tenant compte 
de la reception de I'onde secondaire par 1 'Pigment 
rayonnant (ER) et en mettant en forme en outre : 
10 - une troisiSme matrice d* interaction (F(P')) repr^sentant 
la contribution des sources de 1' obstacle aux points 
predetermines de la surface de 1- element rayonnant (P'l) , 
- et Txne quatridme matrice d' interaction (F'(P')) 
representant la contribution des sources de 1« element 
15 rayonnant aux points predetermines de la surface de 
1 • element rayonnant (P ' i) . 

23. Precede selon I'une des revendi cat ions 19 & 22, dans 
lequel la surface de 1' obstacle correspond t une interface 

20 entre dexix milieux distincts d'une heterostructure . 

24. Precede selon I'une des revendications precedentes, 
dans lequel I'onde principale est line onde sonore et les 
coefficients d6 la matrice d' interaction sont chacun 

25 f one t ion d'un angle d' incidence d'une onde eiementaire 
issue d'une source en ladite region donnee (M) . 

25. Precede selon I'line des revendications 7 et 13, dans 
lequel, pour chaque echantillon de surface, on teste la 

30 valeur d'un produit scalaire entre : 
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- un premier vecteur (r) normal a 1 ' echantillon de surface 
et dirige vers le soramet (P) de 1' hemisphere (fig. 7A) , et 

- un second vecteur {SM) t±r€ entre une source (S) 
associee k cet h^misphSre et ladite region donnge (M) , 

5 en distinguant : 

- le cas ou ce produit scalaire est inf^rieur & un 
seuil pr6d6termin6 et la contribution de cette source 
n'est pas prise en compte, et 

- le cas ou ce produit scalaire est superieur a un 
10 seuil prgdeteirming et la contribution de cette source 

est ef f ectivement prise en compte. 

26, Proc6dg selon I'une des revendications precedentes, 
dans lequel I'onde principale est une onde sonore et, & 
15 I'gtape a), on choisit un nombre total d' echanti lions de 
surface (dSi) sensiblement en fonction d'une longueur 
d'onde de I'onde sonore pour verifier le critSre de 
Rayleigh. 

20 27. Proc6d6 selon I'une des revendications prSc^dentes, 
dans lequel on compare une plurality de valeurs de la 
grandeur physique estimee §l l'€tape d) , obtenues pour une 
pluralite de regions de I'espace, pour selectionner une 
region candidate pour la disposition d'un element 

25 rayonnant destine ii interagir avec 1' obstacle. 

28. Precede selon l'\me des revendications 2 a 27, dans 
lequel I'SlSment rayonnant est un capteur, pour un 
contrSle non destructif, destine t analyser un objet 
30 formant un obstacle de I'onde principale. 
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29. Produit programme d' ordinateur, stocks dans une 
mSmoire d'unit6 centrale ou sur ixn support amovible propre 
Si coop^rer avec un lecteur do cette unit6 centrale, 
caract§ris§ en ce qu'il comporte des instructions pour 
5 mettre en oeuvre le proc§d6 selon l"une des revendi cat ions 
precedentes . 
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